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ABSTRACT 

Qualitative  zesnlts  are  obtained  for  detonation  failnre  due  to  expansion 
loaaes,  The  equations  fo-  the  Mach  number  and  the  reaction  progress 
variable  are  integrated  approximately.  For  certain  loss  ranges  two  det¬ 
onation  speeds  are  obtained.  The  slow  one  is  actually  the  deflagration 
root  that  evolves  beyond  Mach  one.  For  larger  losses  detonation  failure 
oconrs. 


Jilfellnik'JllljlIlM 


SUMMARY  OF  RESEARCH  RESULTS 


To  avoid  repetition,  the  initial  proposal  is  attached  at  the  end 
as  a  part  of  this  Final  Report.  The  basic  equation#  are  Eqs . ( 47) , ( 40) , 
and  (46)  of  the  proposal. 

2(l-m)dm/dz  <1(^  l^q^l/dz 

-  —  •  E  — i— —  —  —  —  — —  — 

m(l+rm)(2+(r-l)tt)  5ixi<ii+Hj+ni2/2 


2e 

( 1+ym) 


(47) 


where  «*’(l/A)dA/dz 


dX^/dz  *  Rj/u 


(40) 


2(y~1)  v  2 

and  c  *=yp/p=  ~~~* - O  >»i«li+ni+1li  ID 

2+(r-l)>»  i 


(46) 


The  crucial  point  is  that  the  right  hand  side  of  eq.(47)  has  to  be  zero  at 

the  sonic  point  whore  the  Mach  number  squared,  M2  =m=l .  This  condition, 

in  conjunction  with  Eq.(40)  will  determine  the  value  of  m  right  behind  the 

2 

initiating  shock  at  z«^0 .  Wo  call  this  m(  z=0)^m1  =  (u]L/c1)  . 


The  above  equations  hove  been  studied  for  realistic  reaction  rates 

1  2 

and  for  different  expansion  lossos.  The  work  of  Tsuge  ,  Wecken  ,  and 

3 

others,  summarized  by  Fickott  end  Davis  ,  has  been  purely  numerical. 
They  Lawo  used  o  thuotin b  to  find  mj,  which  enters  thu 

calculation  as  a  kind  of  eigenvalue.  They  found  that  there  are  two 
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detonation  velocities.  One  increases  and  the  other  one  decreases  as  the 
expansion  loss  increases.  At  a  critical  loss  the  two  roots  coincido. 
Beyond  this  the  two  roots  vanish  indicating  that  there  is  no  steady  deto¬ 
nation  solntion.  This  is  detonation  failnre. 


In  this  work  to  understand  detonation  failure  qualitatively,  I  began 
by  studying  the  phase  plane  of  Eqs.(40)  and  (47).  The  phase  plane  is 
extraordinarily  complex.  The  phase  orbit  begins  at  the  point  z=X=0  where 
m-m^  and,  with  to  be  determined  in  such  a  way  that  the  orbit 
reaches  the  m^l  line  with  a  finite  slope  dm/dX.  This  can  occnr  only  at 
the  intersections  of  the  curves  dm/dX*0  and  m=l .  The  curve  dm/dX*0  is 
not  very  complicated  qualitatively,  even  in  the  cases  where  one  takes 
realistic  temperature  dependent  reaction  rates.  The  problem  is  that  the 
curve  contains  as  a  parameter.  For  different  values  in  the  physical 
range  (Y-l)/2y<m141,  the  curves  look  radically  different;  even  the 
nnmber  of  critical  points  changes  is  one  varies  .  By  this  method  it 
was  not  clear  how  to  go  beyond  the  results  in  Ref. 3  and  categorize  the 
behavior  of  m^  (and  therefore  D)  as  a  function  of  the  loss  parameter 
c •  The  approach  discussed  below  bore  much  more  fruit. 


It  turns  out  thnt  Eq,(47)  can  be  integrated 
integrated  equation  evaluated  at  the  sonic  point  z**zs, 
and  cBcs,  yields  an  exac.  equation  for  m^, 

U-hfOi)2  ^  2 

-  „  (2+(y-l)+2(r-l))  X.q./cj  )A 

(l+yJmj  i 


very  simply.  The 
m=l , X=X8# 


(If) 


where 


esp 


r1  2e 

-  - dz 

J°(l+7«) 
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Vo  get  the  form  of  Eq.(lf),  Eq.(46)  evaluated  at  z=0  has  been  used.  The 
quantities  behind  the  shock  arc  related  to  those  in  front  by  the  usual 
Hugoniot  relations 


2 

S3 


and 


2  +  (y-l)niQ 

2y®0  “  (y-l) 


(2ymg-(y-l))  (2  +(y-1)bq) 
(y+l)2^ 


2 

in  accordance  with  Eq.(lS).  Here  idq^Mq  is  the  Mach  number 

squared  in  the  unreaoted  explosive  in  front  of  the  shock  and  Cg 

is  the  speed  of  sound  there.  The  detonation  velocity  D  is  given  by 
2  2 

D  =®oCq  .  VThen  we  substitute  for  m^  and  Cj  in  Eq.(lf) 
in  terras  of  mg  and  Cq,  we  obtain 

(1  +  y®0^  ^  2 

- (2+(y-l)m0+2(y^l)2 

(1  +  y)^ 

( 2  f ) 

It  may  seem  strange  that  Eq.(2f)  has  the  same  form  as  Eq.(lf),  This  should 
not  surprise  us,  however,  because  the  conservation  equations  are 
integrated  across  the  shock  front  via  the  Hugoniot  relations  for  mg>l, 
whereas  we  simply  put  m^=mg  for  m^l.  The  featuros  of  Eq,(2f)  can 
be  understood  qualitatively.  The  value  of  Dig  is  found  from  tho  inter¬ 
sections  of  the  graphs  of  tho  left  and  right  sidos  of  Eq,(2x)  versus 
mg,  keeping  in  sight  the  Oq  dependence  of  A(c,mg).  For  the 
case  e=0,  Xis=*l,  A(0,mg)=l,  we  obtain  Fig.l.  The  two  roots 
mg_  and  Bg+  are  found.  mg_  represents  a  no-shock  Chapman-Jouguet 
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deflagraion.  m^+  is  the  Chapman-Jouguet  detonation  solution.  As  e 
increases,  ACe.diq)  decreases.  In  the  proposal  the  dependence  of  A 
on  Bq  was  ignored.  Then,  as  e  increases,  the  straight  line  for  e=0 
it  simply  brought  down  with  a  reduced  slope  until  the  two  roots  meet  at 
some  value  m(l  and  then  disappear  below  Mach  one.  The  comments  in  the 
proposal  regarding  Tsuge's  results  were  not  really  valid, then,  because  he 
considered  temperature  dependent  reaction  rates  in  his  numerical 
calculations.  Tho  temperature  dependence  changes  things  beceuse  it  has  an 
important  effect  on  the  dependence  of  A(c,irq).  This  effect 
oan  bo  understood  approximately  by  studying  ACe.m^)  for  a  simple 
case . 


Let  us  take  e  *=  constant  so  it  can  be  brought  out  of  the  integral. 

2 

This  is  the  case  considered  by  Yfecken  (actually,  one  has  to  solve  for 
the  expansion  as  a  free  boundary  problem  of  great  complexity).  The  term 
1  +  ym  in  A  is  a  smooth,  well  behaved  function  of  z.  It  increases 
pretty  much  linearly  from  its  value  1  +  at  z=0  to  its  value 
1  +  y  at  *”*s  as  can  be  seen  by  examining  the  results  in  page  51  of 
Ref. 3.  This  linear  behavior  is  in  marked  contrast  to  the  z  dependence  of 
n,  c.  p,  T,  X,  and  other  relevant  quantities,  which  exhibit  a  character¬ 
istic  step-like  behavior.  Putting  m-1  in  the  integral  for  A  will  not 
affect  the  problem  very  much.  Jn  fact,  putting  m^l  minimizes  the  loss 
effect  of  A.  Therefore  we  will  just  take 


2ez. 


A  -  exp- ( — 


1  +  r 


<  3  f ) 
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The  strong  e»q  dependence  of  A,  due  to  the  temperature  dependent 
reaction  rate,  is  hidden  in  zs.  To  show  this  let  us  take 

dX/dz  =  k  (1  -  X)2  exp(-Ta/T)  . 

where  k  is  a  constant,  for  Eq.(40).  This  is  not  an  entirely  realistic 
reaction  rate  because  k  should  be  a  function  of  T  and  m.  The  T  dependence 
with  the  strongest  effect  on  the  solution  is  contained  in  the  exponential, 
however.  We  write  T  as 

T  =  (M^/Rlp/ p  =  (MM/R)c2/y  , 
or 

2(r-l)q(X+h) 

T  =  (M^/R)  - 

2  +  (y-l)tn 

where  is  the  molecular  weight  of  the  gas,  R  is  the  universal  gas 

constant,  and 


.  c02(  2  +  (y-Duig) 

h  “  (  II.  +  u,  /  2 )  /  q  - -  . 

qfy-1) 

2 

In  the  cases  where  Cq  <<q  wc  can  have  h<<l  if  nig  is  near  one  as 

2 

in  a  weak  detonation.  If  nig  is  large,  however,  h  is  close  to  y  -1 

(y=3  it  a  good  value).  Since  0^X$1  and  m^m^l,  the  strong  detons-- 

tion  behavior  is  well  represented  even  if  we  negluct  the  spatial  varia¬ 


tion  in  T.  Therefore  we  take 


4( y-1 ) q( 1+h) 

T  -  (M  /R) - - - 

(r+i>2 


where  the  minimum  value  y-1/27  has  been,  used  for  m  in  the  denonina- 
tor.  This  overestimate  for  T  makes  the  reaction  rate  as  large  as  it  can  be 
thereby  minimizing  the  loss  effect  dne  to  A.  The  reaction  equation  is 


dX/dz  =  r(l  -  X)2 


with 


( 4  f ) 


r  ■  k  ezp-(b/ ( 1+h) ) 


The  above  expression  can  be  handled  easily  and  has  the  virtue  that  it 
preserves  ti»e  dependence  via  h.  Integrating  Eq.(4f)  with  the  bound¬ 
ary  conditions  X=0  at  z=0,  and  X=Xg  at  z=zs,  we  get 

1-X,  ^  1/ (l+rzg) 

Going  back  to  Eq.(47),  we  impose  the  condition  that  the  right  hand  side  is 
zero  at  X**XS,  «a=l, 

i(l-X8)2  2e 

X,+h  1+7 

(5f ) 

Solving  for  X.  we  get 
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X.  * 


1  - 


2 ( 1 +h) 

1  +  [l+4U+7)(l+h)r/2e]1/2 


(6f) 
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and  then 


*s  “  Xs/r(l-X) 


(7f ) 


The  expression  above  allows  Xs<0  when  h>r ( 1+y) /2r ,  This  is  not 
allowed  and  it  indicates  where  Eq.(Sf)  cannot  be  satisfied.  When  Eqs.(6f) 
and  (7f)  arr  substituted  into  Eqs.(2f)  and  (3f)  one  should  not  take 
seriously  the  portions  where  X$<0.  A  plot  of  Eq.(2f)  is  given  in 
Fig. 2  for  e=0.005  .  In  Fig. 2  the  right  hand  side  has  been  made  zero 
wherever  As<0.  The  two  roots  ihq_=2.74  and  mQ+=10.3  are  found. 

For  this  example  we  have  taken  the  parameters  of  Ref. 3.  page  46  (except 
for  b  which  is  taken  smaller} .  With  these  values  the  solutions  are  well 
into  the  double  detonation  range.  A  natural  d  mensionless  loss  parameter 
6  is  defined  by 

6=(c/k)1/2. 


A  numerical  procedure  was  used  to  find  the  roots.boing  careful  about 
the  k#<0  problem.  The  results  are  presented  in  the  Table  below.  The 
speed  of  sound  in  the  unrcacted  material  is  taken  to  be  Cq*=1367  m/sec. 
at  in  Ref. 3.  The  roots  at  6=0  in  the  Table  arc  the  ideal  CJ  deflagration 
and  detonation  roots.  For  6  less  than  about  0.0025  there  are  still  a 
deflagration  and  a  detonation  root.  A  little  above  that  value  tha  defla¬ 
gration  root  goes  above  Mach  one  to  boc^ne  a  slow  detonation  root.  The 
fast  detonation  root  decreases  while  the  slow  dotonation  root  increases 
with  6.  Near  the  value  6=. 00715  the  two  detonation  speeds  become 
equal  tnd  then  disappear.  For  6=,0072,  there  are  no  roots.  Detonation 
failure  hat  taken  place.  The  variation  of  the  slow  root  is  very  fast  for 
small  values  of  6.  This  is  probably  the  reason  why  Tsuge  could  not 
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follow  the  slow  root  in  his  numerical  calculations.  Thu  values  in  the 
Table  are  plotted  in  Fig. 3.  The  curve  is  very  similar  to  those  obtained  in 
the  numerical  calculations  discussed  in  Ref. 3.  Clearly,  as  the  diameter 
loss  parameter  6  is  reduced,  the  slow  detonation  root  goes  to  Mach  one 
as  conjectured  by  Tsugo^. 


The  above  results  are  based  on  an  approximation  that  overestimates  the 
reaction  rate  by  overestimating  T,  This  givas  a  lower  bound  for  A. 
Fur t hoi  work  is  needed  to  find  an  upper  bound  for  it  and  thus  bound  it 
from  both  sides. 
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TAUI.l-  I 


| 


& 

raO- 

B0  + 

Ds  1  ovr  ( ra/  s  ) 

Dfast(m/s) 

0.0 

0.047 

21  .47 

294.78 

6334.10 

0.000001 

0.052 

21.31 

311.72 

6310.45 

0.00001 

0.067 

20.90 

353 .84 

6249.45 

0.00005 

.0.10 

20.30 

434  .44 

6159.10 

0.0001 

0.14 

19.86 

502.27 

6091.98 

0.0002 

0.19 

19.20 

599.00 

5890.00 

0.0003 

0.25 

18.70 

676.63 

5911.39 

0.0004 

0.30 

18.30 

742 .47 

5847.82 

0.0005 

0.34 

18.0 

797.09 

5799.69 

0.0008 

0.46 

17  .0 

927.14 

5636.29 

0.001 

0.54 

16  .52 

1004.54 

5556.14 

0.002 

0.86 

14  .50 

1268.44 

5205.38 

0.0025 

0.99 

13  .70 

1360.15 

5059.75 

0.0030 

1.35 

12.97 

1588.31 

4923  .10 

0.003  5 

1.71 

12.3 

1787.58 

4794.25 

0.0040 

2.04 

11 .6 

1952.47 

4655 .83 

0.0044 

2.31 

11  .1 

2077  .66 

4554.39 

0.005 

2.74 

10.3 

2262.79 

4387.20 

0,006 

3.60 

8.8 

2593.70 

4055  .18 

0.007 

5.10 

6.9 

3087  .12 

3590.82 

0.0071 

5.40 

6.5 

3176  .62 

3485.18 

0.00715 

5.88 

6  .0 

3314.80 

3351.24 

0.72 

no  root 

no  root 

no  root 

no  root 

L 
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(i) 


-  DESCRIPTION  OF  PROPOSED  RESEARCH 

T-T^rf^  introduction: 

•  •  This  proposal  for  a  mini  grant  is  a  follow-up  to  the  research  the 

"  "  P.I.  performed  at  the  AFATL  at  Eglin  AFB  under  the  1979  USAF-SCEEE  Summer  Faculty 
Research  Program.  That  research  dealt  with  understanding  the  behavior  of  nonideal 
explosives  via  their  mathematical  models.  The  results  of  that  research  are  summa¬ 
rized  in  the  Final  Report  entitled  Detonation  Physics  of  Nonideal  Explosives  With 
Analytical  Results  for  Detonation  Failure  Diameter.  A  copy  of  that  report  is  at¬ 
tached  as  the  Appendix  to  this  proposal. 


The  major  causes  of  nonideal  behavior  discussed  in  the  Appendix  are: 
(i)  endothermic  and  reversible  chemical  reactions,  (ii)  heat  and  a.eral  expansion 
losses,  and  (iii)  transverse  structure  of  the  front.  The  first  two  causes  receive 
some  treatment  in  the  Appendix.  In  the  following  sections  a  continuation  of  the 
research  is  proposed. 


Detonation  Failure: 

The  study  of  cause  (ii)  in  the  Appendix  led  to  the  development 
of  an  analytical  method  that  predicts  detonation  failure  diameters  and  the  exist¬ 
ence  of  a  "slow  detonation"  solution  of  the  steady  state  equations.  A  number  of 
questions  remain  to  ie  clarified.  First,  the  P.I,  believes  that  the  above  solution, 
which  had  been  found  by  other  authors  using  numerical  schemes,  and  which  really 
corresponds  to  a  no-shock  Chapman-Jcuguet  deflagration,  is  an  unstable  solution  un¬ 
der  most  conditions.  In  order  to  demonstrate  this  i~.  is  here  proposed  to  do  a  sta¬ 
bility  analysis.  The  methods  developed  here  will  lo  useful  in  the  work  outlined  in 
the  next  section.  Second,  the  results  in  the  Appendix  rest  on  an  iteration  scheme 
whose  properties  were  not  fully  obtained  quantitatively.  This  was  because  the 
equation  for  the  reaction  progress  variable  cannot  be  integrated  in  closed  form. 
Under  this  grant  it  is  proposed  to  integrate  numerically  Eqs.(48)  and  (61)  of  the 
Appendix  through  several  iterations  to  examine  the  convergence  of  the  iterates 
toward  the  solution  that  is  described  qualitatively  in  the  Appendix.  Since  the 
main  features  of  the  solution  arc  already  understood,  there  is  good  guidance  for  the 
success  of  the  numerical  methods.  The  completion  of  the  above  plan  of  work  would 
be  a  good  start  to  a  truly  2  or  3  dimensional  treatment  of  the  lateral  expansion 
and  the  behavior  of  the  gases  behind  the  sonic  locus. 


(n) 


Transverse  Structure  of  the  Front: 


The  importance  of  cause  (iii)  of  nonideal 
behavior  is  well  established.  A  good  description  of  recent  advances  is  contained 
in  the  bock  by  Fickett  and  Davis^.  The  development  of  the  Transverse  structure 
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It  is  here  proposed  to  use  the  techniques  of  nonlinear  perturbation  and  bi¬ 
furcation  theory  to  extend  the  v/ork  of  Erpenbeck,  cited  in  Ref.  6,  on  the 
formation  of  transverse  front  structures. 


Connection  with  Experiment: 

Prior  tc  my  departure  from  the  HERD  laboratory  at 
Eglin  AFB,  Major  James  W.  Holt,  Jr.,  USAF,  Chief  of  the  Explosives  Branch 
at  the  AFATL,  expressed  an  interest  in  doing  experiments  to  verify  this  theory. 
I  will  reestablish  contact  with  him  to  try  to  arrange  this. 


Facilities  Available: 

The  numerical  calculations  will  be  done  on  the  University 
of  Miami's  UNIVAC  1100/81  computer.  The  budget  of  the  Department  of  Physics 
allows  the  P.I.  to  have  access  to  large  amounts  of  computer  time  without 
charge  to  this  mini -grant. 
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1981  -  82  Est.  salary  $26,620/9  months 

1 00%  time  during  August  3  -  August  15,  1981  (1/2  month)  $1479 

Fringe  Benefits 

14.832  X  faculty  salary  $  718 

Total  of  Salaries,  Wages  and  Benefits  $ 5 5 5 8 


Other  Direct  Costs 
Travel 

Consists  of  transportation  cost  to  scientific 
meetings  and  visits  to  other  Institutions  with 
whom  P.I.  collaborates  on  research  work  per¬ 
taining  to  this  proposal.  $  481 

Publication  Costs  $  200 

Supplies  $  50 

Total  Direct  Costs  $6289 

Indirect  Costs 

(59/  of  modified  total  direct  costs)  $37 1 Q 

TOTAL  $10,000 
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"Entropy  and  the  Approach  to  Equilibrium  of  Infinite  Chains  of  Coupled 
Harmonic  Oscillators",  -Invited  Colloquium,  University  of  South  Florida, 

Tampa,  January  1970. 

2.  "Entropy  Oscillations  and  the  H  Theorem  for  Finite  Segments  of  .Infinite 
Coupled  Harmonic  Oscillator  Chains",  International  Conference  on  Thermo¬ 
dynamics,  Cardiff,  Wales,  April  1970. 

3.  "Exact  Solutions  for  Thermal  Conductances  of  Two  and  Three  Dimensional  Contacts" 
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t  Statement  of  Other  Support  of  Dr.  M.  Huerta 

P.I,  in  AFl  R  Contract  No.  78-3663 

T’  "The  Stability  and  Dynamics  of  Elastic  Structures  and  Fluid  Flows." 

This  contract  supports  the  study  of  imperfect  bifurcations  and  secondary  bi- 
IV  :  atio..-  in  elastic  structures  and  fluid  flows.  It  supports  100C  of  the 
P.I.'s  time  during  May  15,  1980  -  July  6,  1980.  The  budget  for  this  contract 
fol  1  ov/s . 


Revised  Budget  -  AFOSR  78-3663 


Salaries  and  Wages 

Manuel  A.  Huerta 

Principal  Investigator 

1900-31  $24,200/9  months  (est.) 

1 00%  time.  May  15,  1980  -  July  6,  1980  $4705 

Fringe  Benefits  698 

Total  Salaries  and  Fringe  Benefits  $5403 

Travel  Consists  of  transportation  cost  to  scientific 
meetings,  visits  to  and  visits  from  scientists  at  other 
institutions  with  whom  P.I.  actively  collaborates  on 


research  work  pertaining  to  this  proposal  473 

Publication  Costs  350 

Supplies  and  Communications  1 1 0 

Total  Direct  Costs  $6,336 

Indirect  Costs 

( 53 ‘I  of  salaries  and  benefits)  $2,864 

TOTAL  $9,200 
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The  n  reuses  Of  nonldoal  t'chn v < or  In  explosl vos  aro  examined,  An  elabora¬ 
tion  of  thu  2ND  model  for  l-dlmonclonal  steady  detonations  Is  presented  as  a 
reference  frame  from  which  to  examine  the  effects  of  '’■everslbla  chemical 
reactions  endothermic  reactions,  etc,  Mathematical  model*  for  2-dlwonslonal 
detonations  are  studied  and  a  preliminary  analytical  method  of  solution  H 
presented,  The  method  provides  analytical  expressions  for  the  detonation 
failure  diameter  In  terms  of  tho  parameters  of  tho  detonation.  Tho  results 
Obtained  help  to  clarlly  previous  results  obtained  by  other  authors  using 
numerical  methods.  Suggestions  for  further  research  ore  given, 
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I.  INTRODUCTION 

The  most  powerful  high  explosive  formulations  are  usually  composed  of 
TNT,  HMX,  RDX  and  a  few  other  compounds.  These  substances  are  relatively 
expensive  and  their  availability  in  large  quantities  is  not  assured.  To 
head-off  a  possible  shortage  of  high  explosive  materials,  attempts  are  being 
made  by  the  USAF  and  others  to  increase  the  variety  of  compounds  that  can  be 
used  as  high  explosive  components.  Due  to  its  wide  availability,  ammonium 
nitrate  is  one  of  the  most  important  compounds  being  incorporated  into  new 
high  explosive  formulations.  It  is  often  found  that  these  new  formulations 
do  not  perform  as  well  as  one  would  expect  from  their  heats  of  reaction  and 
other  thermodynamic  data.  This  deficient  performance  has  led  to  the  use  of 
the  term  nonideal  explosives. 

There  is  no  hard  definition  of  nonideal  explosives.  The  situation  is 
similar  to  that  for  ideal  vs.  nonideal  gases.  A  nonideal  gas  is  one  whose 
equation  of  state  differs  significantly  from  the  ideal  gas  law  PV=NKT.  No 
real  gas  is  ideal  if  one  looks  hard  enough,  but  under  the  usual  conditions 
of  low  pressure  and  high  temperature,  however,  most  gases  behave  very  closely 
to  ideal  gases.  In  the  case  of  explosives,  ideal  behavior  is  pretty  well 
represented  by  the  Chapman-Jouguet  (CJ)  theory  of  detonations.  In  many  cases, 
where  one  deals  with  explosives  with  fast  reaction  rates  and  large  diameter 
charges,  there  is  good  agreement  between  the  CJ  theory  and  experiment.  In 
many  of  the  new  explosive  formulations  being  developed,  however,  one  finds 
nonideal  behavior.  In  these  no;  ideal  explosives  the  detonation  velocity  and 
pressure  are  found  to  be  significantly  lower  than  predicted  by  the  CJ  theory. 
This  disagreement  between  theory  and  experiment  has  prompted  the  need  for 
further  developments  in  detonation  theory. 

The  Chapman1  -  Jouguot2  theory  was  developed  near  the  turn  of  the  century. 
There  are  many  good  expositions  of  it  sue!)  as  by  Courant  and  Friedrichs3  and 
by  Taylor  and  Tankin'*.  The  CJ  theory  considers  the  reaction  zone  to  be  a  dis¬ 
continuity  in  a  reacting  gas,  very  similar  to  the  usual  theory  of  gas-dynamic 
shocks.  The  fluid-flow  equations  for  ccnrervation  of  mass,  momentum  and 
energy  arc  integrated  across  a  1 -dimensional  plane  detonation  discontinuity. 
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The  above  is  completed  with  the  hypothesis  that  immediately  behind  the  dis¬ 
continuity  the  chemical  reactions  are  completed,  and  the  fluid  velocity  rel¬ 
ative  to  the  discontinuity  is  sonic.  The  CJ  theory  will  not  be  discussed  in 
detail  here  because  more  advanced  developments  are  presented  below.  To  com¬ 
pare  this  theory  with  experiment  one  needs  the  equation  of  state  and  the 
energy  equation  of  the  explosive.  Tluse  equations  are  not  well  known  for 
solid  explosives.  This  is  typical  of  the  present  state  of  the  field,  where 
a  lot  of  important  information  is  just  not  known.  The  situation  for  gases  is 
much  better  than  it  is  for  solids,  and  satisfactory  comparisons  between  theory 
and  experiment  can  be  made. 

The  next  major  step  forward  in  the  understanding  of  detonations  took 
place  during  World  War  II  with  the  development  of  the  Zeldovick-Von  Neumann5- 
Doering  (i.ND)  theory  of  detonations.  A  thorough  discussion  of  this  theory 
as  well  as  many  theoretical  and  experimental  aspects  of  detonations  can  be 
found  in  the  recently  published  book  by  Fickett  and  Davis5.  Von  Neumann's 
rigorous  treatment  of  a  1 -dimensional  steady  detonation  wave  firms  up  and 
clarifies  the  CJ  theory.  His  discussion  of  a  pathological  reaction  also 
introduces  a  way  in  which  nonideal  behavior  can  occur.  His  model  of  a 
detonation  wave  is  that  it  consists  of  a  mechanical  shock  wave  that  initiates 
the  chemical  reaction,  followed  by  a  deflagration  zone  (pressure  drops  away 
from  shock)  where  the  chemical  energy  is  released. 

Present  theories  for  steady  detonation  waves  are  extensions  cf  the  ZND 
theory.  They  include  general  types  of  chemical  reactions  and  deviations 
from  1-dimensional  motion.  These  effects  are  the  main  causes  of  nonideal 
behavior  and  will  be  considered  below.  This  work  is  organized  as  follows. 

Sec  II  states  the  research  objectives.  Sec.  Ill  presents  a  discussion  of  a 
fairly  general  theory  of  steady  detonations.  In  Sec.  IV  the  theory  is 
specialized  to  1-dimensional  detonations  to  highlight  the  nonideal  effects 
due  to  various  types  of  chemical  reactions.  In  Sec  V  the  effects  of  diver¬ 
gent  flow  are  explored.  These  effects  are  quite  complicated  to  calculate 
because  the  flow  is  no  longer  1 -dimensional .  This  section  contains  the 
original  results  on  detonation  failure  obtained  in  this  research  effort. 

Sec.  VI  contains  the  conclusions  and  recommendations. 
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I I .  OBJECTIVES  OF  THE  RESEARCH  EFFORT 

rX-  - 

(1)  To  survey  the  present  state  of  detonation  theory  with  a  viev/  to 
obtain  a  clear  understanding  of  the  causes  of  nonideal  behavior 
in  explosives. 

(2)  To  develop  analytical  methods  for  obtaining  solutions  to  some  of 
the  problems  of  interest.  Specifically,  an  approximation  method 
is  studied  that  can  lead  to  analytical  expressions  for  calculating 
the  detonation  failure  diameter  of  detonating  charges. 

III.  GENERAL  DISCUSSION 

As  the  front  of  a  detonation  wave  passes  a  point,  It  leaves  behind  a  gas 
moving  v/i th  some  velocity.  The  set  of  points  where  the  speed  of  the  moving 

gases  relative  to  the  front  equals  the  local  speed  of  sound  Is  called 
the  sonic  locus.  The  strength  of  th>  detonation  wave  depends  on  the  amount 
of  chemical  energy  that  Is  released  between  the  leading  shock  front  and  the 
sonic  locus.  Any  chemical  energy  released  beyond  the  sonic  locus  Is  lost  for 
the  purpose  of  producing  a  strong  detonation  wave  although  It  may  still  do 
useful  work  against  confining  media.  As  will  be  seen  below,  the  CJ  theory 
predicts  that  in  many  cases  the  sonic  locus  coincides  with  the  point  where 
the  chemical  reaction  is  completed  so  the  entire  energy  released  Is  applied 
to  the  detonation  wave.  Any  deviation  from  this  will  amount  to  nonideal 
behavior.  In  order  to  understand  the  problem  one  should  start  with  a  fairly 
general  model.  A  good  framework  is  given  in  the  papers  of  Kirkwood  and  Wood7 
and  Wood  and  Kii*(-woorJn.  An  even  more  general  starting  point  is  found  In 
Chapters  A  and  5  of  Ref,  G.  The  mathematical  model  uses  the  Euler  equations 
for  a  compressible,  noncJIssIpalivo,  reactive  medium,  undergoing  adiabatic 
motion  into  which  a  steady  detonation  wave  is  propagating.  Let  the  medium 
be  a  cylinder  of  Indefinite  length  hut  finite  radius.  We  Introduce  a  system 
of  cylindrical  coordinates  with  X-axis  along  the  charge  axis,  The  region 
Z  <  0  will  cover  the  unreacted  explosive.  The  shock  front  is  at  Z'Uut  t*0 
and  propagates  toward  The  region  Z  >0  cuvurs  the  zone  of  reaction, 

The  variables  arc  taken  to  be  Independent  of  the  azimuthal  angle  around  the 
Z-axis . 
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The  equations  are: 
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Here  p  Is  the  mass  density;  u  and  <*>  are  the  mass  voloclty  components  In 

the  Z  and  r  directions  respectively;  p  Is  thn  fluid  pressuro;  £  Is  the 

Internal  energy  per  unit  mass;  and  the  are  a  set  of  reaction  progress 

variables  for  the  J  chemical  reactions  taking  place  with  the  reaction  rates 

Rj.  An  equation  of  stato  £  •  t  (p  n ,\ ^ )  Is  assumed  known.  A  form  often 
used  Is  j 

Etf.p;)'  -t--,-  21  Vf;  (o) 

whore  the  q|  are  the  heats  per  unit  moss  released  In  the  reactions.  The 
dot  (.)  notation  means  that  the  connective  derivative  Is  taken 
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We  will  Just  consider  a  steady  detonation  wove  located  at  Zu0,  The  unrcactcd 
fluid  rushes  In  from  the  left  with  a  del onatiun  ipeed  0  (to  he  determined). 

We  also  specialize  the  equations  to  the  1  axis  whore  ruU  and  w**0,  Tho 
result  Is 
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The  Von  Neumann  model  is  to  assume  that  for  Z  <  0  the  flow  quantities 
arc  constant  with  values  pq,  p0,  Dauo>  X^E0.  At  ZB0  there  is  a  mechanical 
shock  with  jumps  to  ,  p^,  but  still  X^O.  For  Z>0  the  reactions  develop 
accompanied  by  changes  in  p,  p,  and  u.  The  Mach  nur, .tiers 
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Where  cQ  is  the  sound  speed  In  front  of  the  shock,  and 


Mi  *  £. 
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Where  Is  the  sound  spped  right  behind  the  shock,  are  related  by 


M*  -  *  +  h^_ 

In  accordance  w.ith  the  usual  shock  relations 
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as  found  in  Landau  and  Lifshitz9.  A  solution  of  tho  shock  conditions  with 
no  discontinuity,  so  P-|  =  PQ,  p-j  =  pq,  and  Mi=M0,  is  also  allowed. 

Eq.  (8)  -  (11 )  can  be  i  ntegrated  with  the  boundary  conditions  at  Z=Q  and  Z=0 


,  Hc')s  ^  '  U,°'T  =  H' 

W0-)=f.  ,  ,  U(o-)-t) 


to  yield  c  =  fu+  2  £.P  i.  - 

and  M.  ‘ID*.  4U|1  '  H<  i  U.1 

where  H  is  the  enthalpy  per  unit  mass, 

H  -  ET  +  P/J 


and 
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,  (21) 

.  (22) 
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in  tho  notation  of  Chapter  5  of  Ref.  6.  The  integrals  I2  and  n  represent 
losses  due  to  the  radially  divergent  flow.  In  one  dimensional  flow  they 
would  be  zero.  Bernoulli's  equation  along  the  axis,  Eq.  (21),  is  not  affected 
by  the  radial  flow.  The  important  equation 
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can  be  obtained  from 
<T: 


Eqs. (0)-(12) .  The  o-  are  given  by 
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The  above  equations  can  describe  the  steady  flow  region. Cq.  (2‘->)  has  a 
singularity  at  u=c,  the  sonic  locus.  Beyond  this  point  the  flow  has  to 
match  the  boundary  conditions  at  Z  large.  Generally  the  flow  beyond  the 
sonic  locus  is  some  sort  of  time  dependent  rarefaction  wave.  Fig.  lisa 
diagram  of  the  situation. 
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Fig.  1.  Main  regions  of  a  detonation  wave  seen  in  the 
frame  of  the  shock  front. 


The  equations  in  this  section  allow  for  a  fairly  realistic  and  mathe¬ 
matically  very  complicated  model.  In  Sec.  IV  we  specialize  the  equations 
to  the  much  simple  1 -dimensional  flow  where  u>=0.  In  Sec.  V  we  consider 
some  aspects  of  2-dimensional  flows. 
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IV.  1-DIMENSIONAL  THEORY 

Putting  w=0  and  3/Dr-O  in  the  equations  of  Sec.  Ill  we  obtain  a 
1-dimensional  theory.  The  results  of  this  section  are  standard  and  are 
presented  to  fix  the  notation  and  nomenclature  that  is  needed  in  Sec.V. 
El  iminating  u  from  Eq.  (19j  and  (20)one  obtains 

H'  '-f.V  ( ‘4  "  4)  -(28) 

In  a  p  vs.  1/p  plot,  the  graph  of  Eq.  (28)  is  a  straight  line  called  the 
Rayleigh  line.  Eliminating  u  and  D  from  Eq.  (21)  one  obtains 

■  >/f).  <2S> 


which  is  a  family  (as  one  varies  )  of  curves,  in  a  p  vs.  Vp  plot,  known 
as  Hugoniot  curves.  As  we  move  in  Z  behind  the  front, the  A.  are  related 
because  from  Eq.  (12)  we  have 
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Thus  in  effect  a  single  A  parametrizes  the  Hugoniot  family.  So  far  D  is  not 
determined  so  the  slope  of  the  Rayleigh  line  (-p2D2  )  is  unknown.  Yet  the 
intersection  of  the  Rayleigh  line  with  the  Hugoniot  curves  will  determine 
the  solution.  The  answer  is  obtained  by  noting  that  Eq.  (25)  gives  a 
singularity  for'3u/az  at  the  sonic  locus  unless  we  have 

■23^  =  o  •  on 

(we  are  in  1 -dimension,  so  w=0).  Eq.  (31)  must  be  imposed  at  the  sonic  locus 
to  avoid  the  singularity  for  the  reasons  given  in  Chapter  5  of  Ref.  6. 

Eq.  (31)  determines  the  value  of  the  A.  and  thus  picks  a  particular 
Hugoniot  curve,  at  the  sonic  locus.  Noting  that  l)-c  at  the  sonic  locus,  and 
using  Cq.  (19),  we  have  that  the  slope  of  the  Rayleigh  line  is 
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(32) 


We  also  have  that  the  slope  of  the  Hugoniot  curves  is 
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Combining  Eqs.  (27),  (32),  and  (33),  we  get  for  the  slope  of  the  Hugoniot 
curves 


M  -  -  fV 
Vf>  I 


Thus,  the  slope  of  the  Rayleigh  line  is  the  same  as  the  sIodp  of  the  Huqoniot 
curve  (v/hich  is  determined  by  Eq.  (31)  )  at  the  sonic  locus.  For  the  case 
where  there  is  only  one  reaction  Eq.  (31  )  becomes  R=0  or  x  =  l  at  the  sonic 
locus.  This  means  that  the  one  reaction  is  complete  at  the  sonic  locus. 

This  justifies  the  CJ  hypothesis  v/hich  determines  the  Rayleigh  line  (and  thus  D) 
by  putting  it  tangent  to  the  completed  reaction  Hugoniot.  For  the  equation 
of  state  of  Eq.  (5),  the  above  method  yields  that  0  satisfies 

D*  -  f  >  v  +  cl  =  o  .  (35) 

L  C?-  J 

Eq.  (35)  has  the  tv;o  roots 
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The  two  roots  give  two  Rayleigh  lines  called  R+  and  R-.  R+  represents  a  CJ 
-detonation  while  D-  represents  a  CJ  deflagration.  These  are  depicted  in 
Fig.  2.  For  the  CJ-  solution  the  solution  goes  smoothly  along  R-  from  to 
the  values  at  CJ-  attained  at  the  sonic  locus.  No  shock  takes  place  for 
this  deflagration  solution.  For  the  detonation,  the  Von  Neumann  model  is 
that  the  solution  goes  along  the  no  reaction  Hugoniot  from  0  to  N, 
p=:p^and  p=o^.  (with  Pi  and  p^  determined  from  Eqs.  (15)  and  (17)  using 
Eg  (36)  for  MQ)  and  then  down  from  N  to  CJ+  along  R+,  A  detonation  solution  going 
directly  from 


->% 


Fig.  2.  The  p~  V  p  plane.  CJ+  gives  the  pressure  at  the  sonic  locusfor  a 
detonation  .CO -  gives  the  pressure  for  a  deflagration 


0  to  CJ+  along  R+  with  no  shock  is  also  allowed.  Both  solutions  have  the 
same  final  state  CJ+  at  the  sonic  locus  and  the  same  D.  The  no  shock  solu¬ 
tion  does  not  usually  occur  because  of  problems  with  initiating  the  reaction. 
The  Von  Neumann  model  in  fact  says  that  the  0-*-CJ<  detonation  is  actually  a 
mechanical  shock  to  pj  followed  by  a  deflagration  from  N  to  CJ+.  Due  to  the 
way  that  the  equations  integrate  through  a  shock,  the  point  CJ+  is  the  some 
with  or  without  a  Von  Neumann  spike. 


The  CJ  hypothesis  is  that  the  sonic  locus  is  determined  by  tangercy  of 
the  Rayleigh  line  to  the  completed  reaction  Hugoniot.  If  there  are  several 
Irreversible  exothermic,  reactions  this  hypothesis  is  correct  (for 
1 -dimensional ,  steady  flows).  If  some  of  the  reactions  are  reversible,  however, 
or  if  some  of  the  a  j  in  ay  be  negative  as  in  mole  decrement  reactions  ,  then 
Eq.  (31  )  docs  not  imply  that  the  reactions  are  completed  at  the  sonic  locus. 
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Thus  some  of  the  chemical  energy  is  released  beyond  it,  This  late  energy  is 
not  effective  in  the  dr  loin  Li  on  .md  nonidr.il  behavior  is  found.  The  analysis 
of  precisely  what  happens  involves  a  study  of  the  phase  plane  of  Eq.  (25). 

Many  complicated  results  are  possible  depending  on  the  details  of  the o ^  and 
the  R^ .  A  fairly  thorough  discussion  of  the  possibilities  can  bo  found  In 
chapter  5  of  Ref.  6,  and  in  Wood  and  Salsburg11  . 

This  brief  discussion  of  steady,  1 -dimensional  detonations  docs  not 
cover  the  entire  subject.  In  this  research  effort  we  have  mado  no  con¬ 
tribution  to  this  part  of  the  field.  There  are  good  analytical  methods, 
that  supplemented  by  computer  calculations,  can  predict  the  noriidoal  per¬ 
formance  of  steady,  1 -dimensional  detonations.  An  Input  to  those  methods, 
however,  must  contain  a  realistic  description  of  the  reaction  kinetics.  This 
Is  presently  not  known  for  most  gaseous  explosives,  let  alone  solid  ones.  In 
the  case  of  solid  explosives,  beyond  the  lack  of  information  concerning  tho 
equation  of  state  and  the  reaction  kinetics,  there  arc  also  the  problems 
associated  with  Inhomogeneous  mixtures.  Here  the  effects  duo  to  grain  burning 
and  diffusion  of  reactants  are  not  fully  worked  out. 


V.  DIVERGENT  FLOW  AND  DETONATION  TAILUIIE 


When  u  Is  not  negligible,  Egs.  (g)-(l2)  need  to  be  complemented  by  a 


way  to  compute  as  a  function  of  1.  This  depends  on  tho  shape  of  tho 
exploding  charge  and  on  the  properti es  of  the  confining  media,  In  general, 
to  compute  w  one  must  solve  a  difficult  free  boundary  problem  with  matching 
at  tiie  boundary  to  account  ior  the  properties  o (  tho  different  media.  Even 
if  u  were  known,  however ,  there  are  still  some  interesting  difficulties  in 
integrating  Eqs .  (d)-(|2)-  These  equation',  should  predict  the  observed 
phenomenon  of  detonation  failure  when  the  diameter  of  the  exploding  cylinder 
is  small.  We  v/ill  cast  those  equation:,  in  the  form  for  flow  in  a  nozzle  by 
letting 
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where  A  is  the  area  of  the  nozzle.  The  equal  inns  become 
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Several  authors  hove  token  the  above  equation-,  with  realistic  rate  terns  and 
applied  computer  methods  to  obtain  a  variety  of  results.  Some  of  this  work 
1*  tuning  Hied  In  Chaplet  £  of  Kef.  L.  Wockon"  takes 

1  A  A  „  p  _  constant  ,  (42) 

A  fT7 

with  one  Arrhenius  reaction  to  discuss  the  structure  of  the  reaction  zone. 

He  is  not  successful  In  obtaining  detonation  failure  at  some  critical  c  because 
of  computational  difficulties.  Tsugn  et  a  1 . 1 7  take 

di. ,  -i  [  yi-.-.? _ ]Vi .  (43) 

a  v  ^  \_.(  \  v6Vt  )l\  -P‘J  J 

where  Is  the  value  of  the  thickness  u  of  an  explosive  slab  at  Is  Q,<j  Is 
given  In  Cq.  (1C)  , 

A  -  <L  ,  Pi  li-A  ,  «nd  It  -  41l  •  (44) 
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where  PQ0  is  the  density  of  the  confining  medium.  Eq.  (43)  is  obtained  by 

treating  the  lateral  expansion  as  a  blunt  body  about  which  the  confining 

Inert  flows.  Tsuge  et  al.  take  a  very  detailed  description  of  the  reaction 

kinetics  for  H2  +  1_  02->-ll?0  and  solve  the  problem  numerically.  Their  ccm- 
2 

parisons  with  experiment  for  exploding  slabs  are  fairly  good.  Their  numeri¬ 
cal  calculations  lead  to  plots  as  shown  in  Fig.  3.  For  c-|<oc  there  is  no 

detonation  solution,  so  detonation  failure  is  found  in  this  model.  For  c  >o 

1  c 

they  obtain  two  solution  branches.  The  upper  branch  is  like  the  usual 

1-dimensional  detonation  M  of  Eq.  (36).  They  did  not  explore  the  lower 

branch  fully  because  of  computational  ai fficul ties .  They  speculate  that  as 

°1‘>  ~  the  lower  branch  has  M^l  .  It  will  be  shown  here  that  this  is 

incorrect.  The  limit  of  the  lower  branch  as  o  *►  ®  is  actually  M  -*■  M 

1  o  o- 

as  in  Eq.  (36 ) , 


/or 


Rather  than  attacking  this  problem  numerically  from  the  beginning  we  wil 
try  to  develop  the  solution  analytically  as  for  as  possible-  Vie  use  the 
equation  of  state  given  in  Eq.  (6).  T linn  the  o..  defined  in  Eq.  £6)  (not 
to  be  confused  with  the  u  and  above  which  were  lengths)  are  given  by 


FT,  =  U-_0  ,  (40) 

O' 

where  , 

^  *|2  --  T\u  +  U  v  u ,l/Ll  .  (46) 
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We  rewrite  Eq.  (39)  as 
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where  /Y>\=  ^  =  W/ci  .  (47) 


l  4-xvn 


v  art 
_  > 


A 


Eq.  (47)  is  to  be  solved  so  that  dm/dz  remains  finite  at  m=l •  A  necessary 


m=1 .  Also  at  Z=0,  \=0  and  m=m^  with  to  be  determined 


condition  for  this  is  that  the  right  hand  side  of  Eq.  (47)  be  zero  when 

These  are 

unusual  boundary  conditions  to  impose  on  a  first  order  ordinary  differential 

equation.  Eq.  (47)  does  not  satisfy  the  usual  Lipschitz  condition  needed  for 

the  existence  of  a  solution  with  a  desired  value  of  m  at  a  specific  Z  .  But 

the  unusual  boundary  conditions  are  sufficiently  "free"  that  a  solution  will 

exist.  We  do  not  fix  the  Z  where  m=l  occurs,  the  sonic  point  Z,  but  we  let 

the  equation  tell  us  so  as  to  keep  dm/dz  finite.  Then  we  find  the  value 

•  Since  the  problem  presented  by  Eq.  (47)  is  not  a  standard  one  it  is 

unclear  that  a  Picard  iteration  method  will  converge  for  it.  We  have  found 

useful  results  with  this  method,  so  we  will  use  it  with  the  awareness  that 

t  h 

it  may  need  modification  due  to  convergence  difficulties.  The  nL  iteration 
of  m ,  mn  will  be  found  from 


-1 
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where  n  =  2,3,.;.. 
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-  ~  ,  (48) 
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with  the  right  hand  side  of  Eq.  (AO)  equal  to  zero  when  mn  -  1 .  We  also 
need  to  use  Eq.  (AO).  The  first  iterate  m1  is  found  by  putting  the  dA/dz 


term  equal  to  zero.  The  equation  foi 


in1 


can  be  integrated  to  get 


(18) 


,  (49) 


rV[2-vCt-ow'J  Ll+Sf  ^<1 


cr 
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where  >yi't  -  W(£ -o)  ,  (50) 

and  where  we  have  used  Eq.  (46)  evaluated  at  Z--0  ,  with  c-j  =  c(Z=0). 
Imoosing  the  right  hand  side  condition  for  the  m  equation  we  get 

.  4. 

±  ~6  at  TTi  =  l  .  (51) 

If  all  the  qi  are  positive,  and  the  are  as  given  in  Eq.  (5),  we  need 
that  for  the  first  iteration, at  the  sonic  point  Z‘  , 

7Y\(t\)~  1  ,  .  (52) 

For  the  simple  this  implies  Z*  -  »  .  We  go  to  Eq.  (49)  and 

put  *.j  =1  when  m1  =1  to  find  m|  .  T'-.e  result  is  that  mj  obeys 

the  equation 

o  4-  3(v-0  .  (53) 

(tv*)'4-  ’rn!  c?  ' 


Manipulating  this  equation  yields 


(54) 


(19) 


.Eq.  (54)  is  almost  the  same  as  Eq.  (35).  The  difference  is  that  in  Eq.  (35) 


the 


we  have  cQ  and  0=  MQcQ  whereas  appears  in  Eq.  (54).  If  we  consider 

.  ,  -  m1  ,  and  we  have  the  same  result  as 

01  0 


case  with  no  shock,  then  r  cQ,  mj 
In  Eq.  (35).  The  solutions  obtained  would  refer  to  a  no-shock  detonation  that 
goes  along  0-CJ+  and  a  no-shock  deflagration  that  goes  along  0-CJ-.  The 
second  of  these  is  meaningful,  but  not  the  first.  If  we  consider  the  case  with 
a  shock  we  put 

C*  V  I5* 

?•  (y  v Yn0 


(55) 


Into  Eq.  (54)  to  get,  after  using  Eq.  (15), 


(56) 


This  is  exactly  the  same  as  Eq.  (35),  just  as  it  should  be  because  the  R+ 
line  in  Fig.  2  is  the  same  whether  we  consider  a  shock  solution  along 
0-N-CJ+  ,  or  a  no  shock  one  along  0-CJ,  .  The  two  solutions  of  Eq.  (55)  are 
interpreted  as  a  Von  Neumann  detonation  along  0-N-CJ^  ,  and  an 

unphysical  shocked  deflagration  along  0-S-CJ-  .  Out  of  the  four  solutions 

obtained  we  only  desire  the  shocked  detonation  0-ll-CJ+>  and  the  no-shock 
deflagration  0-CJ-  .  Their  detonation  velocities  are  given  by  the  roots 

D+  =  c0MQ+and  D_  =  C0M0_  of  Ec)»  (35).  Using  Eq.  (15)  we  get  the  mj  + 

that  goes  with  M  while  the  m,1  that  goes  with  fl  Is 

m]_  =  •  ■  Knowing  mj  +  and  m|_  we  go  to 

Eq.  (49)  and  we  get  mj  and  m*  as  functions  of  X^z).  Finally  from 
Eq.  (  40)  vie  can  integrate  to  get  Xl(z)  .  This  first  iterate  reproduces 
the  results  of  Sec.  IV  because  dA/dz  y/us  neglected.  But  we  begin  to  sec 
the  source  of  the  two  roots  in  the  work  ofTsutjc  et  al.1?  When  we  compute 
the  next  iteration  the  results  of  the  expansion  will  come  in.  But  as  the 
diameter  goes  to  "  the  expansion  becomes  negligible  and  the  results  reduce 
to  and  •  To  compute  nr  we  go  to  Eq.  (4P)wlth  n»2. 

He  integrate  it  with  the  boundary  condition 


(20) 


Yr\\ 
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The  result  Is 
Jn 


0  +  •*  T a'-Y  r>\\  [l  M  V  -0»n]  tf%  ^  % 

J  V  / 

(50) 


m|  Is  determined  by  Imposing  the  bounded  dm?dz  condition  to  tho  right 
hand  side  of  tq.  (48)  with  n-2,  At  the  sonic  point  zj  w#  have 
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Eq.  (40)  was  used  to  obtain  Eq.  (GO).  The  lattnr  equation  says  that  ^  f  ]  at 
the  sonic  point  so  tho  roactlons  ire  not  completed.  Some  of  the  energy 
will  bo  roloasod  lata  leading  to  nonldeal  behavior,  In  principle  tq.  (59) 
determines  l\  in  terms  of  mj  .  This  can  be  soon  as  follows.  We 
haw  that 


M,  4  (WiV  6  4=  (£}1  I  ,  (GO) 

~Z  L  l,‘ 

•-  ft 


50  that  cln  l.q.  (GP  )  Is  known  In  {ntm%  of  m*  end  Cj  .Also  from  Cq.  I  40) 

T *i  C 


(ni) 
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In  Eq.  (61)  m2  is  known  in  terms  of  a2  from  Eq. (58)  •  (We  assume  that 
A  dA/dz  is  explicitly  known  as  a  function  of  Z  ).  Therefore  one  can 
integrate  Eq.  \61)  to  get  ,\2  as  a  function  of  Z  and  m^ 

When  this  function  A2(  Z,  m2  )  is  substituted  into  Eq.  (59)  one 

can  find  Z2  as  a  function  of  m2  as  claimed.  Now  we  go  back  to 
Eq.(58)  and  demand  that 

^(Vs)=  ViU  =  X  .  (62) 


and  get 

0+V)  v\* 


In  Eq.  (63)  Z2  ana 
equation  for  m|  . 

A  =■  ^ 


are  functions  of 


.  So  Eq.  (63)  is  an 
Tfie  last  terms  is  the  expansion  dependent  term.  Let 
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Here  e  represents  an  inverse  diameter  that  goes  to  zero  as  the  diameter  goes 
to  ”  •  Also  a(c)  -*■  0  as  e  0  or  as  the  expansion  is  negligible. 

Eq.  (63)  is  rewritten  as 


+»=o  ■ (65) 


This  equation  reduces  to  Eq.  (54)  as  a->  0  and  A2  -*•  1  as  e  0. 

Thus  there  is  a  clear  connection  between  the  two  branches  to  be  obtained  from 
Eq.  (65)  and  thel-dimensional  results.  The  two  roots  of  Eq.  (65)  are  obtained 
'  by  the  quadratic  formula  as 


(66) 


(22) 


e  »’ 


For  small  e  these  two  roots  give  the  mj  and  mi  of  before 

:^;.5&which  were  connected  with  the  mi+  and  m*_  of  Eq.(35).  As 
,  e  and  A  grow  beyond  a  critical  value  the  two  roots  become  complex,  so 
4  there  is  no  real  solution  and  detonation  failure  takes  place.  To 
simplify  things  let  us  ignore  the  variation  of  A2  with  e  • 

Then  a  plot  of  the  roots  versus  A  appears  as  in  Fig.  4.  This  is  qualitatively 
the  same  as  the  curve  ^  ^ 


obtained  numerically  by  Tsuge  et  al.For  a  <  Ac  (for  the  diameter  greater 
than  some  critical  value)  there  are  two  branches.  The  upper  branch  is 
connected  as  A-*o  to  the  1-dimensional  shocked  detonation  0-N-CJ+  of  Fig.  2. 
The  lower  branch  is  connected  as  A**0  to  the  1 -dimensional  no-shock 
deflagration  0-CJ_  .Now  one  sees  the  error  in  the  work  Tsuge  et  al.  They 
could  compute  things  like  Z2  and  ^  numerically  but  with  so  much 
difficulty  that  the  results  were  hard  to  interpret. 

The  qualitative  features  of  the  solution  are  now  clear.  To  get  quanti¬ 
tative  results  one  needs  to  integrate  Eq.  (61).  This  would  have  to  be  done 
numerically.  It  is  clear  that  successive  steps  in  the  iteration  are  not 
very  different  from  the  steps  to  calculate  m2  .  No  theoretical  problems 
with  convergence  should  arise  although  this  should  be  checked  numerically 
in  some  specific  examples.  We  propose  to  do  this  in  future  work.  The 
critical  value  A  at  which  detonation  failure  occurs  is  given  by 
I  -i-  Of-0  ^  [i  v  Z'KV  -  C  Y  .  (67) 

It  would  be  interesting  to  explore  ways  in  which  this  behavior  could 
be  checked  with  experiment.  It  should  be  noted  that  this  theory 

predicts  that  the  speed  of  deflagrations  increases  as  the  diameter  decreases. 
At  some  critical  diameter  the  detonation  and  deflagration  velocities  are  the 
same  and  for  smaller  diameters  there  is  no  solution. 


(23) 


VI  CONCLUSIONS  AND  RECOMMENDATIONS 

The  causes  of  nonideal  behavior  can  be  summarized  in  three  categories.  ^ 

(i)  Incompletely  effective  use  of  the  chemical  energy  for  the 

detonation  wave  due  to  the  nature  of  the  chemical  reactions. 

If  there  are  reversible  or  endothermic  reactions,  the  effects  an 
especially  important.  For  steady  1 -dimensional  flows  the 
physical  tfieory  is  prepared  to  handle  the  problem.  The  knowledg 
of  the  burning  rates  of  the  relevant  chemical  reactions,  however 
needs  to  be  improved. 

(ii)  Loss  of  strength  due  to  the  lateral  expansion  of  the  explosive. 
Here  the  state  of  the  physical  theory  is  not  so  good.  In  Sec.  V 
we  presented  a  calculation  for  detonation  failure.  The  radial 
flow  was  not  calculated  but  was  assumed  known.  The  analytical 
solution  presented  should  be  checked  with  numerical  calculations 
A  further  step  would  be  to  compare  these  results  with  experiment 
The  author  would  like  to  be  able  to  pursue  these  further 
developments. 

(iii)  Loss  of  strength  due  to  transverse  structure  of  the  front  and 
non-steady  phenomena  in  the  detonation.  These  effects  are 
interesting  per  se  and  their  effects  on  nonideal  behavior  are 
judged  to  be  important6.  These  wave  instability  problems  were 
not  addressed  in  this  research  effort  because  the  author  felt 
he  needed  a  good  understanding  of  the  simpler  aspects  of  noni¬ 
deal  behavior  before  addressing  then.  There  is  much  room  for 
research  into  these  effects  using  the  techniques  of  nonlinear 
stability  analysis  and  bifurcation  theory. 


